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ABSTRACT 

The gamma-ray space telescopes AGILE and Fermi detected short and bright synchrotron gamma-ray flares 
at photon energies above 100 MeV in the Crab Nebula. This discovery suggests that electron-positron pairs 
in the nebula are accelerated to PeV energies in a milliGauss magnetic field, which is difficult to explain with 
classical models of particle acceleration and pulsar wind nebulae. We investigate whether particle acceleration 
in a magnetic reconnection layer can account for the puzzling properties of the flares. We numerically integrate 
relativistic test-particle orbits in the vicinity of the layer, including the radiation reaction force, and using 
analytical expressions for the large-scale electromagnetic fields. As they get accelerated by the reconnection 
electric field, the particles are focused deep inside the current layer where the magnetic field is small. The 
electrons suffer less from synchrotron losses and are accelerated to extremely high energies. Population studies 
show that, at the end of the layer, the particle distribution piles up at the maximum energy given by the electric 
potential drop and is focused into a thin fan beam. Applying this model to the Crab Nebula, we find that the 
emerging synchrotron emission spectrum peaks above 100 MeV and is close to the spectral shape of a single 
electron. The flare inverse Compton emission is negligible and no detectable emission is expected at other 
wavelengths. This mechanism provides a plausible explanation for the gamma-ray flares in the Crab Nebula 
and could be at work in other astrophysical objects such as relativistic jets in active galactic nuclei. 
Subject headings: Acceleration of particles — Magnetic reconnection — Radiation mechanisms: non-thermal 
— ISM: individual (Crab Nebula) — Gamma rays: stars 
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1. INTRODUCTION 

The conversion of magnetic to particle kinetic energy is a 
long standing problem in astrophysics, e.g. in relativistic jets 
and pulsar winds (the "er" problem). Magnetic reconnection 
is the main mechanism known to dissipate magnetic energy 
into thermal and non-ther mal energy of the plasma (see e.g. 
IZweibel & Yama da 2009 for a recent review). A precise and 
accurate understanding of this process is still lacking, in par- 
ticular for relativistic reconnection (i.e. if the Alfven speed 
approaches the speed of light, V\ ~ c), but current Particle- 
In-Cell (PIC) simulations have become a powerful tool to 
probe the microphysics of reconnection. Several of these 
PIC simulations were carried out in the context of relativistic 
electron-positron pair plasmas as found, for i nstance, in pul- 
sar w i nds (see e.g. iZenitani & Hosh ino 2001; larosche k et al.l 
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Liu et all \20l ll ISironi & Spitkovskvl 1201 lb . These studies 



demonstrated that magnetic reconnection is able to channel 
a sizeable fraction of the available magnetic energy into a 
population of non-thermal ultrarelativistic particles via the 
induced large-scale electric field. 

The gamma-ray space telescopes AGILE and Fermi recently 
discovered bright day- long gamma-ray fl a res above 100 MeV 
in the Crab Nebula jAbdo et al.l 1201 It iTavani et al l 1201 It 
iBalbo et al.l l20TTt IStriani et al.l 1201 U iBuehler et al.l 1201 lF . 
This high-energy emission is thought to be synchrotron radia- 
tion by relativistic electron-positron pairs in the nebula. How- 
ever, these flares exhibit several puzzling prope rties that chal- 
lenge classical models of pulsar wind nebulae dRees & Gunnl 
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l!974tlKennel & Coronitil 1984b and particle acceleration. One 
of the most intriguing features is the emission of synchrotron 
photons above the classical limit = (9m e c 2 / '4ap)(E / B ±) ss 
160 (E/B±) MeV (where m e is the mass of the electron, a? 
is the fine structure constant, E is the electric field and B± is 
the magnetic field perpendicular to the particle's motion) im- 
posed by the balance between the accelerating electric force 
and the dampi ng radiation-reaction force due to synchro tron 
energy losses dGuilbert et al.|[T983t Ide Jager et al.l 1 1996b . as- 
suming that E < B±. These observations support the presence 
of PeV particles in the nebula (corresponding to Lorentz fac- 
tor 7 e ~ 10 9 ), the highest energy particles ever asso ciated with 
a specific astrophysical source dAbdo et all 1201 lb . Another 
challenging aspect of the flares is that most of the observed 
gamma-ray emission above 100 MeV should originate from a 
tiny part of the nebula (cffl are ~ 3 x 10 15 cm <C 0.1 pc) if the 

emitting region is causa lly connected. 

In a previous study (lUzdenskv et al.l 1201 lb . we proposed 
that the gamma-ray flares in the Crab Nebula could be at- 
tributed to extreme particle acceleration in the vicinity of a re- 
connection layer present in the nebula. In this scenario, high- 
energy particles move across the reversing magnetic field and 
follow the relativistic analog of Speiser's trajectories dSpeiserl 
11965b . We found that the particles' orbits are naturally fo- 
cused into the reconnection layer where the reconnecting 
magnetic field is small, in particular smaller than the recon- 
nection electric field. Because of this, particles suffer weaker 
synchrotron losses and can be accelerated to extremely high 
energies deep inside the layer (Kirk 2004). Motivated by the 
results of our first study, here we examine in detail the process 
of particle acceleration in the electromagnetic field configu- 
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FIG. 1. — Top: Diagram of the magnetic field line configuration (isocon- 
tour of the magnetic flux function ip) at the reconnection site in the .ry-plane 
containing the reconnecting (B x ) and the reconnected (B y ) components of the 
magnetic field. The reconnection layer is delimited by the grey region be- 
tween y = ±<5. Bottom: Same as in the top panel, but where the reconnection 
layer is broken up into a chain of magnetic islands. 

ration likely to be found in a reconnecting system. We use 
relativistic test-particle simulations, i.e. the interaction be- 
tween particles and the back-reaction of particles on the fields 
are neglected. Importantly, in the equation of motion we in- 
clude the radiation reaction force induced by synchrotron and 
inverse Compton losses. The only other sim ilar studies that 
also included the radiation reaction f orce are IS chopper et alJ 
(fl998l fT999h . and lNodes etaflll200l using resistive magne- 
tohydro dynamics (MHD) simulation s; and the PIC simula- 
tions by iJaroschek & Hoshinol d2009t) . In this paper, we fo- 
cus specifically on the motion of the highest energy particles 
(with a Larmor radius greater than the thickness of the recon- 
nection layer) present in the plasma and discuss the possible 
application of this process to the gamma-ray flares in the Crab 
Nebula. 

This article is organized as follows. In Section|2j the main 
assumptions and approximations adopted for the reconnection 
layer in a relativistic pair plasma are presented. Section [3] is 
dedicated to the equations of motion of relativistic particles in 
the reconnection layer, with a particular emphasis on the ex- 
pression for the radiation reaction force. Section@]focuses on 
the motion of a single particle and analyzes the effect of each 
large-scale electromagnetic field component, and synchrotron 
and Compton losses. We study the distribution of a population 
of electrons at the end of the reconnection layer in Section|5] 
and apply the model to the Crab Nebula in Section [6] The 
final part (Section|7]) summarizes the results of this study. 

2. THE RELATIVISTIC RECONNECTION LAYER 

2.1. Emphasis on the highest energy particles 

In the present study, we are interested in the dynamics of 
the highest energy particles present in the reconnection re- 
gion, namely, the particles with a Larmor radius much larger 
than the thickness of the layer. These particles are assumed 
to be preaccelerated by some unspecified mechanism. In the 
context of the Crab Nebula gamma-ray flares, the emitting 



particles have a typical Lorentz factor 7 e > 10 9 and a gyro- 
radius p e ~ 10 15 cm (in a milli-Gauss field), of order the size 
of the emitting region L ~ c?fl are ~ 3 x 10 15 cm that we asso- 
ciate here with the length of the reconnection layer L z . The 
emphasis on such particles is unique, because one does not 
usually face the situation where p e ~ L z in other astrophysical 
environments, e.g. in the solar corona. Because of their large 
Larmor radii, the particles of interest here feel only the large- 
scale field structure as described below. Small-scale (e.g. tur- 
bulent) structures in the current layer are probably important 
for acceleration of lower energy particles, but do not strongly 
affect the motion of the highest-energy particles. This justi- 
fies the assumptions of using only the large-scale fields in this 
study, and significantly simplifies the problem. 

2.2. The large-scale electromagnetic fields 

We consider a simple geometry of a steady state and lam- 
inar magnetic reconnection region in a collisionless plasma 
of relativistic pairs (Fig. [T). The large-scale electromagnetic 
field consists of the reconnecting magnetic field B x in the x- 
direction, the reconnection electric field E z in the z-direction, 
and the reconnected magnetic field B y in the y-direction. In 
the general case, there is also a component of the magnetic 
field along the electric field B z , called the guide field, and 
the associated in-plane induced electric field. The guide 
field B z has a major effect in suppressing relativistic drift- 
kink instabilities in the reconnection lay er and thus may pro- 
mote non-thermal particle acceleration (IZenitani & Hoshinol 
120081) . Because there is no Hall effect in pair plasmas, the 
quadrupole out-of-plane Hall magnetic field and the associ- 
ated electrostatic in-plane electric field components (impor- 
tant in electron-ion reconnection) are absent. 

The current layer has a slab-like geometry with a character- 
istic thickness S in the y-direction across which the magnetic 
field lines reverse and reconnect. In collisionless pair recon- 
nection, the lower limit on the thickness of the layer is set by 
the relativistic electron skin depth d e = (7 av m e c 2 /47r« e e 2 ) 1 / 2 
(where « e is the electron density and 7 av is the characteris- 
tic Lorentz factor of the particles in the plasma) if there is no 
strong guide field. Otherwise, 5 is given by the Larmor radius 
of the bulk electrons 8 = 7 av po, where po = m e c 2 /eBo. These 
two scales are comparable if the guide field strength is similar 
to (or smaller than) the reconnecting magnetic field and if the 
upstream plasma pressure is not large compared with the mag- 
netic pressure. Note that if the layer is very long (compared 
with its thickness 5), it is subject to instabilities (e.g. tearing) 
that may lead to its effective broadening. In this case, and as- 
suming that the plasma is incompressible, the effective thick- 
ness o f the layer is at most 5 = (3 rec L x (e.g. IShibata & Tanumal 
120011) . where f3 rec < 1 is the dimensionless reconnection rate 
and L x is the size of the layer in the x-direction. 

The reconnecting magnetic field B x is roughly uniform out- 
side the reconnection layer (\y\ ^> 6) and equals ±Bq. Inside 
the layer, the amplitude of B x decreases and by symmetry van- 
ishes at y = 0. The field is commonly modeled by a Harris 
profile, 

flx(y)=Botanh(|). (1) 

The electric field induced by magnetic reconnection is 
roughly uniform in the reconnection region in a steady state, 
and is a finite fraction of the reconnecting magnetic field: 
E z = firecBo- The reconnection rate (3 KC is poorly constrained 
in the context of collisionless relativistic pair plasmas but a 
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FIG. 2. — The blue arrows represent the electric field vector in the xv'-plane 
induced by the motion of the plasma through the reconnection layer in the 
presence of a guide field. 



few s imulations estimate /3 lec ~ 0.1 (e.g. IZenitani & Hoshinol 
120081) . We will adopt this value in this paper. 

In the presence of the guide field, the motion of the plasma 
in the reconnection layer ge nerates an electric field in the xy- 
plane dSchopper et al.ll 19991) . such that E = -V x B z /c. Here, 
V is the velocity of the plasma in the xy-plane which we model 
as (assuming that the Alfven speed Va « c) 



y * =c UJ cosh u 

V y =-c/3 rec tanh(|j . 



(2) 
(3) 



These components of the electric field E x and E y are included 
in this study (see Fig.|2]for a visual representation of the elec- 
tric field structure in the xy-plane). Any bulk motion of the 
plasma in the z-direction is neglected, V z = 0. 

The structure of the reconnected magnetic field B y is much 
more uncertain than that of the other large-scale fields. The 
reconnection region can be viewed as a long Sweet-Parker- 
like layer of length L z ;§> S and width along the x-direction 
2L X , which is generally similar to L z . In the layer (|x| < L x ), 
the large-scale reconnected magnetic field behaves approxi- 
mately as 



B y (x) = (3 lec B Q [ — 



(4) 



reaching P kc Bq at the end of the layer, \x\ = L x . Outside the 
reconnection layer (\x\ > L x ), B y increases rapidly with |x| to 
reach Bq over a distance comparable to L x . 

2.3. Tearing and kink instabilities in the layer 

This simple picture may be a reasonable approximation 
for the average large-scale field structure in the reconnection 
region, but it may not capture its small and intermediate 
scale structu re. In the Sweet-Parker picture (ISweetl IT958b 
lParkerl[T957l) . the expected reconnection rate would be very 
small, much sm aller than 0.1. Petschek-like reconnection 
(Petschek 1964) would be much faster but this configuration 
is unlikely to appear in pai r plasmas because there is no 
Hall effect (see, however, Bessho & Bhattacharjeel 120051 
120071) . Analytical calculations and numerical simulations 
show instead that a thin Sweet-Parker-like layer is unstable 
to secondary tearing modes and breaks up into a highly 
dynamical chain of magnetic islands (or plasmoids) con- 
nected by small current layers (see e.g.. | Loureiro et al1l2007b 
ISamtanev et al.ll2Q09t iDaughton et al.ll2009ir" In this regime, 
magnetic reconnection is fast (see e.g.. IShibata & Tanumal 
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be 0.1 (IZenitani & Hoshinol 120081) . It is important to note 
that each magnetic island contributes to B y , but the net re- 
connected magnetic flux within an island is zero because the 
magnetic field lines are closed loops (Fig. Q] bottom panel). 
The effective large-scale reconnecting magnetic field would 
result from the semi-open reconnected field lines generated 
betwe en the current layers and the islands dUzdenskv et al] 
l2O!0l) . 

A simple way to quantify the effect of magnetic islands on 
the propagation of particles is to add a static perturbation if; to 
the undisturbed magnetic flux function ipo (this quantity cor- 
responds to the z-component of the magnetic vector potential 
A z ) of the form 



ifi = eBoScos 



«i S 7T I — 



(5) 



where «i s is the number of islands and e is the amplitude of 
the perturbation normalized to BqS, both quantities consid- 
ered here as free parameters. Only the y-component of the 
magnetic field is perturbed, and equals 



dtp 

B y =-— = en is irB \ — 



(6) 



The overall smoothed growth of the reconnected magnetic 
field would not be very different from a linear variation across 
the layer as described by Eq. @. 

A thin current layer can also bend along the z- 
direction due to the rela t ivistic drift kin k instability (e.g. 
IZenitani & Hoshinol 120051 : iLiu et al.ll201 ll) . unless there is a 
strong guide field. We describe the effect of the kink insta- 
bility as a stationary sinusoidal translation of the current layer 
in the y-direction as a function of z, so that the reconnecting 
magnetic field is given by 



B x (x,y) = Botanh 



y n ■ {2nn z z 



(7) 



where C y is the amplitude of the deformation of the layer in 
the y-direction and n z is the number of oscillations along the 
layer. This simple prescription is a rather good representation 
of the layer in the linear stage of the kink instability. The other 
electromagnetic fields are unchanged. 

3. RADIATION REACTION FORCE AND EQUATIONS OF MOTION 

This section gives the complete expressions for the equa- 
tions of motion of a single electron moving in an arbitrary 
elec trom agnetic fi eld and s ubject to radiative damping (Sec- 
tion l3.41 >. Sections [3 . 1 113 .21 focus on general and app roxim ate 
expressions of the radiation reaction force. Section [33] pro- 
vides the expression of the inverse Compton drag force in the 
Thomson regime. 

3.1. General expression for the radiation reaction force 

Independently of the underlying acceleration mechanism, 
a relativistic particle of charge e and Lorentz factor 7 e loses 
energy at a r ate given by th e Larmor formula (in covariant 
form, see e.g. I.Tacksonl[T975l) 



(8) 
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where = (7 e ,7 e v/c) is the four- velocity, and the interval 
ds = cdt /7e (t is the time in the laboratory frame). This expres- 
sion can be reformulated in terms of the components of the 
three-acceleration a parallel (an) and perpendicular (a \ ) to 
the di rection of motion of the particle as dRvbicki & Lightmanl 
fT979h 

2e 2 



3c 3 7e 



(9) 



In return, radiative energy losses affect the dynamics of the 
particle through the radiative back-reaction four-force g M = 
(7eg-v/c 2 ,7 e g/c) given by dJacksonll 1975b 



2e 2 d 2 u> 
3c ds 2 



ad 



(10) 



The first term in this equation is known as the Schott term 
and the second term describes the recoil of the particle due 
to radiative losses. In the ultrarelativistic limit (7 e ;§> 1, the 
regime we are interested in here), the Schott term is negli- 
gible and the spatial com ponent of g 11 can be expressed as 
(lLandau & Lifshitd[l97lt) 



rad 



(11) 



and acts as a friction force opposite to the direction of motion 
of the particle. If the radiative losses term vanishes, T^d = 0, 
the Schott term determines the strength of the radiation reac- 
tion force. But then, the radiation reaction force is negligible 
compared with the electric force in the equation of motion 
(see below). 

3.2. Radiation reaction force in an electromagnetic field 

If one now considers the specific situation of a relativistic 
electron of mass m e moving in a given electromagnetic field 
(E,B), the equation of motion is (Lorentz-A braham-Dirac 
equation in the laboratory frame, lJackson|[l975l) 



m e c- 



du^ 
ds~ 



c 



(12) 



where F^ v is the external electromagnetic field-strength 
tensor. In the framew ork of classical electrodynamics, 
lLandau & Lifshitzl (11971b derived an approximate expression 
for g in terms of the external electromagnetic fields using an 
iterative procedure. First, the radiation reaction four-force is 
neglected in the equation of motion, g^ = and only the ac- 
celeration due to the Lorentz force contributes to the radiative 
energy losses. Substituting the four-acceleration term du^/ds 
by -{e/m^F^u,, in Eq. ® yields 



<r> 2 2 

>>rad= 3 r e c 7e 



E+- 



vxB 



v-E 

c 



(13) 



where r e = e 2 /m e c 2 is the classical radius of the electron. The 
term proportional to B 2 in Eq. ([T3l is the classical synchrotron 
energy loss related to the perpendicular acceleration by the 
magnetic field (a± « eB±/j e m e , B± is the component of the 
magnetic field perpendicular to the particle's direction of mo- 
tion). This is the dominant term in most cases. The other 
terms are associated with the linear acceleration by the elec- 
tric field and the perpendicular acceleration by the electric 
field and the magnetic field. These terms are included in our 
calculations although we find that they do not lead to signifi- 
cant changes. 



In the first order approximation for j e 3> 1, the radiation 
reaction force is 



g = 



ad 



-V. 



(14) 



The higher order terms in the expression of the force can be 
obtained by successive iterations (by injecting into Eq. (O 
the Lorentz force and the first order radiation reaction force 
term). In the second order (in r )~ 2 ), the radiation reaction force 
contributes to the linear acceleration term in the Larmor for- 
mula (Eq.[9]l such that 7 e fl||. ra d = |g|/7 e w e («± is unchanged). 
Hence, the first order formula in Eq. (fPfl) is valid as long as 
a± > 7efl|i,rad, implying that 



7 e fi 



(15) 



i.e., if the magnetic field in the rest frame of the particle 
does not exceed the classical critical magnetic field strength 
B c = m 2 ,c 4 /e 3 ~ 6 x 10 15 G (above which the formal Larmor 
radius of the particle in the rest frame p e = m e c 2 /eB becomes 
smaller than r e ). In reality, quantum effects become important 
at a much smaller field Bqed = m 2 ^ /he = otpB c rj 4.4 x 10 13 G 
(at which p e is comparable to the reduced Compton wave- 
length A e = h/m e c) and so the condition in Eq. ( fT3b should be 
reformulated as b = 7 e B /Bqed -C 1 . In the quantum regime 
(b ^> 1), the electron emits in discrete steps photons wi th an 
energy comparable to the electron's energy (lShenlll972h and 
the radiation reaction force should be modeled as a stochastic 
process where the electron underg oes discrete and subs tantial 
energy losses A^ e /j e w 1 (see e.g. iDuclous et al]l20l"R) . This 
situation would require a full quantum electrodynamics treat- 
ment which is not undertaken in this article. 

It is worthwhile to note that the radiation reaction three- 
force can be comparable or even exceed the Lorentz three- 
force Fl without violating the condition b <C 1 . Indeed, with 
|g| ~ r eleB 2 and |Fl| ~ eB, the ratio of the two forces scales 



as 



|f l | 



B c 



= a F 7e&- 



(16) 



Hence, if |g| ~ |F L | and b <C 1, then the electron has to be ul- 
trarelativistic with a Lorentz factor 7 e > ap 1 3> 1. In the con- 
text of the gamma-ray flares in the Crab Nebula, where PeV 
electrons (7 e ~ 10 9 ) evolve in a milli-Gauss magnetic field 
(see Section|6]l, b ~ 10~ 7 and we are well within the range of 
validity of Eq. ( TBi i even when the particles reach the radiation 
reaction limit regime where |g| ~ |Fl|. 

3.3. Inverse Compton drag force in the Thomson regime 

If the electron is bathed in an external radiation field, it 
loses energy via inverse Compton scattering. In the Thomson 
regime (e /m e c 2 <S 1, where e is the target photon energy in 
the rest frame of the electron), the energy loss in each scatter- 
ing event is small, A7 e /7 e ^ 1, and the Compton drag force 
can be described as a continuous force. If the photon field 
is isotropic in the laboratory fr ame, then the Compton pow er 
radiated by a single electron is (iBlumenthal & Gouldlll970l) 



"Pic = -^ot:c%U+ 



(17) 



where ctt = (8/3)irr 2 is the Thomson cross section and £4 is 
the energy density of the radiation field (in erg cm" 3 ). In the 
rest frame of the particle, the target photons are focused in a 



Extreme particle acceleration in magnetic reconnection layers 



5 



cone of semi-aperture angle ~ 1 /j e in the direction opposite 
to the motion of the particle. Hence, the Compton drag force 
can be formally written as a friction force such that 



(18) 



The continuous drag force approximation is not valid in the 
Klein-Nishina regime (e /m e c 2 3> 1) in which the electron 
loses a significant fraction of its energy in a single interaction 
(as for synchrotron losses in the quantum regime). 



Particle' s orbit 




3.4. Equations of motion 

Summarizing the discussion in the preceding sections, the 
equations describing the time evolution of an ultrarelativistic 
electron are (in the laboratory frame, and for b <C 1) 



dt \ c 
df 



Vic 



(19) 



(20) 



In these equations p = j e m e \ is the momentum and £ e = 
j e m e c 2 is the total energy of the electron. In the following, 
we will refer to Eqs. (fT9ll20b as the equations of motion. 

4. MOTION OF A SINGLE ELECTRON IN THE MAGNETIC 
RECONNECTION CONFIGURATION 

We solve the equations of motion including the radiation 
reaction in the vicinity of the reconnection layer described 
in Section [2] We examine in detail the effect of each com- 
ponent of the magnetic field on the dynamics of a single 
particle, starting with the reconnecting component B x only. 
We first review the m ain soluti ons of the analytic model of 
lUzdenskv et al.l (1201 ll) (Section l4~Tl i. Analytical solutio ns ar e 
then compared with the numerical solutions (Section I4.21 i. 
We then investigate the effects of the guide field components 
(Section |4~3l l, the reconnected field and t he p resence of mag- 
netic islands in the current she et (Section |4~4"i i, and kinking of 
the current sheet (Section [4.5b on the motion of the par ticle . 
The effect of Compton scattering is discussed in Section l4~6l 

4.1. Analytic solutions in the reconnecting fields only 

lUzdenskv et al.l (1201 ll) studied the evolution of a single 
electron, assuming that the guide field and the reconnected 
magnetic field are negligible (and without magnetic islands), 
and neglecting the inverse Compton drag force. In this model, 
the electron is initially injected at the origin of the axes with 
a Lorentz factor 70 3> 1 and a velocity contained in the yz- 
plane crossing the z-direction with an angle Qq < ir/2 such 
that v y (z = 0) « c sin Qq (see Fig. [3). If the Larmor radius of 
the electron is much greater than the thickness of the layer, 
the trajectory of the electron is composed of portions of cy- 
clotron orbits stretched along the z-direction by the reconnec- 
tion electric field. The reversing reconnecting magnetic field 
always deflects the particle towards the midplane y = 0. As 
the electron gains energy, the Larmor radius of the particle 
increases but the maximum deviation of the orbit from the 
midplane y max decreases. 

It is possible to derive simple analytical solutions to the 
equations of motion Eqs. (fT9]|20l in three special regimes: 

(1) Most of the particle's trajectory is well outside the 
reconnection layer (y max 3> S) and its energy is well below 



FlG . 3 . — Geometrical quantities defining the orbit of a single electron in the 
vicinity of the reconnection layer (relativistic Speiser orbit) in the yz-plane. 

the radiation reaction limit (ecEo ^> V Ta d), so that radiative 
losses can be ignored. Taking the y-component of Eq. (T% 
and writing d/df = v z d/dz and defining ojq = eBo/m e c, we 
obtain (over one half-cycle, where y > 0) 



dz ^ 



-wo- 



(21) 



Integrating this equation yields 

7e/3 y =-z+7o/? y ,o, (22) 

where z = z/po, Py = v y /c and By : o = /3 y (0) « sin O - Consider- 
ing, for simplicity, the case 9q <C 1, so that /3 Z = v z /c ~ 1 and 
8 y ~ dy/dz (with y = y/po), Eq. d22l i becomes 



dy 

7e-7= =-Z+J(A). 

dz 



(23) 



The evolution of the Lorentz factor of the particle follows (af- 
ter integration of Eq. ( f20b ) 

7e® = fteeZ + 70- (24) 

Combining Eq. d23l and Eq. d24"i i, the final integration over z 
gives the orbit of the electron over one half cycle 



y(z) = ~ + (jf) (1 +/3 re c0o)ln ( 1 + ^ 

ftec VAec/ V 70 



(25) 



If the amount of energy gained by the particle over one half 
cycle is small such that A7 e /7o = /3 rec z/'Jo -C 1, then the tra- 
jectory of the electron follows (after an expansion of Eq. d25l l 
to third order in /3 rec f/7o) 



z 2 1 8 z 3 

y® « 6oz- — ( 1 + + - ^ . 

270 3 7q 



(26) 



At the end of the half cycle, the particle crosses the z-axis at a 
distance approximately equal to 



L 2 
Zl W 2706*0+ -70Prec#0- 



(27) 



The variation of the Lorentz factor and the midplane crossing 
angle over one half cycle are then 



A7e = 7e (Zl ) - 70 ~ 2^ rec 7or? 
A I ^0 I = I ^0 (Z 1 ) I - ^0 « - ^ /? re c ^5 - 



(28) 
(29) 



The evolution of the crossing angle 9q over many cycles 
(z > Zi) is governed by d|6\)|/dz « A|6> |/zi and hence 
after integration one finds that \8q\ oc (7 e )~ 2//3 oc (z)~ 2//3 . A 
similar derivation gives that the maximum distance from 
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In this 



the midplane decreases as y max oc (7 e ) -1 ^ 3 oc (f) -1 ^ 3 . 
regime, the energy of the particle increases and the particle's 
orbit shrinks steadily towards the midplane. 

(2) The particle reaches the radiation-reaction limit 
(ecEo = Had) while well outside the layer (y max 3> S). 
This means that the energy gained by the electric acceleration 
is completely radiated away by the electron in each cycle. 
The Lorentz factor of the electron over each cycle is constant 
and is determined by the balance between the acceleration 
rate and the synchrotron energy loss (neglecting the energy 
loss by the electric field in Eq. (fOll), i.e. 



7e = Trad : 



3c^ rec 
2r e cj 



1/2 



(30) 



This is the standard radiation reaction limit, resulting in the 
synchrotron photon energy limit at ~ 160 (Eq/Bq) MeV intro- 
duced in Section [T] However, in spite of the constant energy 
of the electron, the orbit continues to shrink towards the mid- 
plane because of the radiation reaction force. Using the same 
notation as in regime (1) and replacing 7 e with the expression 
for 7 ra d, the ^-component of the equation of motion becomes 



7rad 

After integration (at t ■■ 
1), we have 



dz 

0, /? y (0) 



l-/3rec/3y 



(31) 



i sin#o « 9q, assuming 9q <C 



1 + /3rec/3y (z) = ( 1 + AeA) exp 



7rad 



(32) 



Expanding this solution to the second order in /3 re cz/7rad <C 1, 
and integrating over z, the equation of the trajectory becomes 



7 2 1 B 

■^(l+/3rec^0)+7 

270 6 



Tad 



The trajectory intersects the y = plane at a distance 



Zl 



27,-ad^O- ^7rad/3rec6'o- 



(33) 



(34) 



A 1 a 



(35) 



Over one half cycle, the midplane crossing angle varies by 

2 

The long term evolution of 9q over many cycles follows 
d#o/dz ~ A 1 9o \/z\, hence the midplane crossing angle 
decreases exponentially as #0 oc exp (- /3 lec z / 3j m d) ■ 

(3) In the third case, one considers an electron whose 
orbit is entirely confined deep inside the layer (y max <C 6) 
where the magnetic field varies approximately linearly 
Bx(y) ~ (y/S)Bo. The maximum magnetic field strength felt 
by the particle is B x (y max ) = (y mlix /S)Bo. Assuming that the 
electron remains in the radiation reaction limit within each 
cycle, the maximum energy reached increases as the magnetic 
field decreases such that 7 e = (<5/y ma x)7rad- The joint evolution 
°f #0' ymax an d 7e m this regime can be roughly estimated if 
one assumes that this situation is similar to the first regime 
described above where Bq is replaced by (y max / 8)Bq. Then, 



one finds 8q oc 7, 



-3/2 



t [light days] 
0.04 




It is important to note that a particle will not always 



FIG. 4. — Example of a numerically calculated relativistic Speiser orbit 
(blue solid line in panel (a) of an electron initially injected in the yz-plane 
at the origin of the axes with a Lorentz factor 7i„j = 5 X 10 s and 8 m j = 90°. 
Outside the layer, the reconnecting magnetic field is Bo = 5 mG. The electric 
field is uniform with /3 1CC = 0.1 and the length of the layer hi = 0.1 light- 
day. The red dotted lines delimit the thickness of the reconnection layer, y = 
±6 = ±10 6 p - The bottom panel (b) compares the evolution of the radiative 
losses (red line) with the almost constant electric acceleration rate (blue line). 
The green dashed line is the radiative loss averaged over each cycle. The 
reconnected B y and the guide B z fields are neglected. 

undergo the three regimes in a sequence. This will depend 
on the initial conditions imposed on the electron and on the 
properties of the reconnection layer. 

4.2. Comparison with numerical solutions 

In order to track the particle's orbit accurately over many 
cycles across the reconnection layer, a high-order numeri- 
cal scheme is required. We use the e xplicit Runge -Kutta- 
Verner procedure of the eighth order ( IVernerf[l978l) . Fig- 
ureHJshows an example of a numerically integrated orbit of an 
electron evolving in a uniform electric field E z = 0.1 Bq with 
Bq = 5 mG and the reconnecting magnetic field B x reversing 
over a layer of finite thickness S = 3.4 x 10 1 1 cm (correspond- 
ing to S = p e (7e = 10 6 )). The reconnected B y and the guide B z 
fields are neglected at this stage in order to compare with the 
above analytical solutions. As expected, we find that the orbit 
shrinks steadily toward the midplane. 

Figure [5] shows the long-term evolution of the particle's 
orbit parameters r y t , and 6q. The initial parameters of the 
electron were chosen specifically to obtain all three regimes 
described in the previous section in a sequence. At first 
(regime 1, z < 0.75 light-day, 1 light-day « 3 x 10 15 cm), 
most of the electron's orbit is well outside the layer (y max 3> <5, 
8 = 3.4 x 10 11 cm) and the electron undergoes almost pure 
electric acceleration (7 e ^ 7 ra d)- The midplane crossing an- 

—2/3 

gle almost follows the power-law 9q oc 7 C . The deviations 
from the analytic solution are entirely attributed to the effect 
of the radiation reaction force, not completely negligible here. 
Then, the electron enters the radiation-reaction limit regime 
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FIG. 5. — Panel (a): Evolution of the electron's Lorentz factor ■y e (red solid 
line) and the midplane-crossing angle 6q (blue dotted-dashed line) along the 
orbit inside the reconnection layer of length L z = 10 light days = 3 X 10 16 cm 
and thickness 6 = 10 6 po = 3.4 X 10" cm. Panel (b): Evolution of the electric 
acceleration rate (blue line) and the cycle-averaged radiative losses (green 
line). Panel (c): Relation between 8q and 7 e calculated along th e orb it. The 
green dashed lines give the analytical solutions derived in Section |4~Tl (power- 
laws of index —2/3 and —3/2). In all panels, the dashed vertical line indicates 
the place where the orbit of the electron lies completely inside the current 
layer (y m ax < The gray region show the period where the electron is in 
the radiation reaction limit regime (-y e = 7 ra d). So = 5 mG, /3 r ec = 0.1 and 
the reconnected and the guide fields are neglected as well as Compton drag 
(Z4=0). 

(regime 2, 0.75 < z/(l light day) < 2.75, gray shaded area) 
where the Lorentz factor saturates at 7 rac j ps 4.3 x 10 8 and 
Oq decreases exponentially (linearly on a linear-logarithmic 
scale), in agreement with the analytical model. After some 
time, the orbit shrinks inside the layer (regime 3, y max < S, 
Z> 1.2 X 10 16 cm). Radiative losses no longer counterbalance 
the electric acceleration exactly over each cycle and the ratio 
of these two quantities tends to about ~ 2/3. The electron is 
continuously focused and accelerated in the z-direction to en- 
ergies above the standard radiation reaction limit associated 
with the B field (7 e > 7 raa ), and 9 decreases asymptotically 



with energy as 0q oc 7, 



-3/2 



4.3. The guide field 

In addition to B x and E z , we study in this section the effect 
of a uniform finite guide field B z and the induced in-plane 
electric fields E x and E y on the motion of the particle. The 



reconnected field B y and magnetic islands are still neglected 
at this stage. Figure [6] presents a few numerically integrated 
orbits projected onto the yz-plane where only the guide field 
strength changes from B z = to B z = 0.5 Bo and B z = Bq. We 
find that the presence of a strong guide field (B z > Bo) sup- 
presses the focusing of the particle toward the reconnection 
layer, even if the particle is injected inside the layer at a small 
angle, 6q <C 1. For B z = Bo, the particle's orbit is still tied to 
the layer (the particle does not escape in the y-direction) but 
the particle remains well outside of the layer (j max ;§> 8) most 
of the time where its energy is limited by the radiation reac- 
tion force. For a moderate guide field strength (B z = 0.5 Bo), 
the particle focuses toward the midplane as in the case with- 
out guide field, but at a slower rate (see the middle panel in 
Figure |6j. 

Examination of the three-dimensional orbit for B z = Bo 
shows that the particle spirals along the guide field in the xy- 
plane and drifts slightly in the x-direction (Fig. [7). This drift 
in the x-direction as well as the defocusing of the particle's 
orbit in the yz-plane in a strong guide field is mainly due to 
the in-plane electric field E x . Indeed, E x accelerates the parti- 
cle coherently in the x-direction as the particle makes one turn 
along the guide field in the xy-plane, resulting in an increase 
of the Larmor radius of the particle. In contrast, E y results 
only in a small net energy gain because it first decelerates and 
then accelerates the particle in the y-direction by almost the 
same amount during one cycle along B z (depending on where 
the particle is in the xy-plane, see Fig. [2j». To summarize this 
section, we find that a strong guide field B z > Bo, does not 
allow the particles to converge into the layer. As a result, 
particles cannot be accelerated above the standard radiation 
reaction limited energy 7 ra( j. 

4.4. The reconnected field and magnetic islands 

We next investigate the effect of the reconnected magnetic 
field By on the motion of the particles in the vicinity of the 
layer. We will first assume that B y follows the linear law given 
by Eq. (0]i and magnetic islands are still neglected at this stage. 
To single out the effect of B y only, there is no guide field and 
no Compton drag in this subsection. All fields are fully in- 
cluded later in the population studies (Section|5}. 

In Figure [8] we calculate the Lorentz factor of the particle 
7 0u t and the angle of deviation of the orbit 6> out from the z- 
axis when the particle reaches the boundaries of the layer (i.e. 
when z = L z or |x| = L x ), as a function of the initial injection 
distance xq (between xq = and xo = L x , yo = and zq = 0). 
The initial particle velocity is along the y-axis with a Lorentz 
factor 70 = 10 7 . The layer is L z = 2L X = 4 light-day-long and 
S = 3.4 x 10 11 cm thick, with Bq = 5 mG. As expected, we 
find that the reconnected magnetic field B y deflects the par- 
ticle in the x-direction away from the z-axis. The deviation 
angle 6> out increases quasi-linearly with the injection distance 
xo because of the linear increase of B y with x. The break 
at Xofrk/Lx ~ 0.18 is geometrical in origin because it distin- 
guishes whether the particle escapes the layer at z = z ou t = L z 
or at x = x out = L x . Below xo.bk, the particle exits the layer in 
the z-direction first (x out < L x ), and above x~bk the particle is de- 
flected out of the layer in the x-direction first (z ou t < L z ), This 
explains why the particle is accelerated to the highest energy 
if xo < Xo bk because it goes through the total electric potential 
drop available, e(3 KC BoL z « 2 PeV (see Fig. [8] bottom panel). 
In addition, we observe that the presence of the reconnected 
field B y does not affect the focusing mechanism of the parti- 
cle's orbit toward the midplane, but instead of having a beam 



8 



Cerutti, Uzdensky, & Begelman 




10 

o 

-10 



I B z /B o =0.5 



B =0 




I B Z /B D =1 




FIG. 6. — Examples of electron orbits projected in the yz-plane for a finite guide field B z =0 (top), 0.5 Bo (middle), and Bo (bottom). The particle is initially 
injected with 70 = 10 7 at the origin along the v-axis (80 = 90°). The layer is L z = 3 X 10 15 cm ( 1 -light day) long, of thickness 8 = 3.4 X 10" cm, with Bo = 5 mG, 
fte C = 0.1,B y = 0andW*=0. 




FIG. 7. — Three-dimensional electron orbit corresponding to the case where 
B z = Bo in Fig. [5] The grey lines show the projections of the orbit in the yz- 
and .ry-planes. 

of particles focused in the z-direction, the trajectories spread 
over a wide range of angles |0 out | < 90° in the xz-plane at the 
end of the reconnection layer (fan beam of particles). 

As discussed in Section [2] a thin current layer is unsta- 
ble to tearing modes, leading to the formation of a dynam- 
ical chain of secondary magnetic islands (where fields lines 
are closed loops: see Fig. [1] bottom panel, for an illustra- 
tion). In order to quantify the impact of magnetic islands in 
the layer on the dynamics of the particles, we add to the lin- 
early growing reconnected field B y a sinusoidal perturbation 

By prescribed in Eq. (O (see the discussion in Section|2). Is- 
lands are static and identical flux ropes, of length L z in the 
2-direction and of width W x = 2L X /«; S in the ^-direction. In 
the v-direction, the width of each island is about W y ~ eS if 

e > 1, and W y ~ (2e)^ 2 5 if e < 1 . 

Calculations show that in some cases, the particle is trapped 
within one island and its orbit stays aligned along the z- 
direction. In the xy-plane, the particle's orbit is confined and 
focuses toward the center of the island (O-point). Figure [9] 
shows that the three-dimensional trajectory of the particle be- 




0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 8. — Lorentz factor (bottom panel) and deviation angle (defined with 
respect to the z-axis, top panel) of the electron at the end of the reconnection 
layer as a function of the initial position of the particle inside the layer xq < 
L x . Initially, the electron is injected along the y-axis (6q = 90° ) with a Lorentz 
factor 70 = 10 7 , inside a layer of length L z = 2L X = 1.2 X 10 16 cm (4-light 
days) and thickness S = 3.4 X 10" cm, with Bo = 5 mG and ft cc = 0.1. The 
vertical dotted line separates the region where the particle escapes the layer 
in the z-direction first (xq/L x < 0.18) from the region where it escapes in 
the .v-direction first (*o/L x > 0.18). The horizontal red dashed line marks 
7max = eprecBgLz/ (m^c 2 ). There are no magnetic islands and B z = in this 
calculation. 

comes quite complex and it is rapidly trapped by one island. 
For this calculation, we considered the same particle (70 = 
10 7 ) as in Figure[8]with xq = for a L z = 1 light-day-long layer, 
S = 3.4 x 10 11 cm and Bo = 5 mG. The layer contains n; s = 500 
islands and e = 1.4 so that B y , max = 0.5 Bo- This magnetic 
confinement by the islands operates if the large scale recon- 
nected field (linearly increasing with x) is locally smaller than 
the amplitude of the island field B y < B y . max = eni s irBo(S/L x ), 
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FIG. 9. — Top panel: Example of a three-dimensional orbit of a particle 
trapped in a magnetic island present in the reconnection layer. The layer is 
L z = 1 light-day-long, S = 3.4 X 10" cm with Bo =5 mG and contains a static 
chain of «j 8 = 500 islands with e = 1.4. The particle is injected at the origin 
along the y-axis with 70 = 10 7 . There is no guide field. The grey lines show 
the projections of the orbit in the vz- and .vy-planes. Bottom panel: Projection 
of the orbit in the .vy-plane overlaid on the magnetic field lines of two islands 
in the layer. 

i.e. if the particle is injected at a distance xq < eTrn^S / (3 lec . 
Otherwise, the particle feels only the large scale magnetic de- 
flections by By. If the amplitude B y , max exceeds the maximum 
large scale reconnected field By. max = PnxBo, then we find that 
the particle is systematically confined within an island for any 
value of xo, and regardless of the size of the island with re- 
spect to the Larmor radius of the particle. We conclude that 
the presence of magnetic islands in the current sheet does not 
hamper the acceleration of the high-energy particles of inter- 
est here; they even help in keeping particles aligned along the 
electric field. Note however that this simple model is static, 
i.e. it does not take into account any electric fields associated 
with the contraction or the merging of islands, which could 
affect particle acceleration in a s imilar way as for smaller par- 
ticle energies dDrake et al.ll2006l) . 

4.5. Kink instability 

We now consider a sinusoidal deformation of the current 
layer in the y-direction as a function of z, using the expres- 
sion in Eq. (jT) for the reconnecting field to mimic the ef- 
fect of a kink instability (during the linear or early nonlinear 
development). The calculation of particle trajectories shows 
that the deformation of the layer, even with high amplitudes 
(C y > 1), does not change at all the focusing mechanism. The 



final energy reached by the particle at the end of the layer is 
also unchanged. Figure [10] is an example of a particle trajec- 
tory initially injected at the origin along the y-direction with 
70 = 10 7 in a layer with C y = 1 and n z = 5. In this calcula- 
tion, the layer is L z = 1 light-day long and 6 = 3.4 x 10 11 cm 
thick with Bq = 5 mG and /? rec = 0.1. There are no magnetic is- 
lands, guide field or external photon fields for simplicity. The 
layer guides the particle trajectory. The kink affects the en- 
ergy of the particle and the focusing mechanism only if the 
wavelength of the kink mode considered X z = L z /n z becomes 
shorter than the distance traveled by the particle alo ng th e z- 
direction during one half Speiser cycle z\ (see Sect. I4.lt . In 
this case (A z <C z\), the particle always feels a non-zero aver- 
age perpendicular magnetic field during each cycle. Its energy 
becomes limited by the radiation reaction force and the parti- 
cle never reaches the inside of the layer. 

4.6. Inverse Compton cooling 

The inverse Compton drag force strongly affects the trajec- 
tory of the electron if the radiation energy density is of order 
the magnetic energy density U+ > Ub = B^/Sir. To quantify 
the effect of inverse Compton drag, we compute the deviation 
angle out and the Lorentz factor 7 out of the particle at the end 
of the layer as a function of the ratio U+/Ub- FigurefTTIshows 
the result of this calculation where the particle is injected at 
different locations xo in the midplane (yo = and zo = 0). The 
layer has the same properties as in Section l4~4l with B z = and 
no magnetic islands. For simplicity, the background of target 
photons is homogeneous and isotropic in the reconnection re- 
gion. We find that inverse Compton cooling increases the de- 
flection of the particle's trajectory by the reconnected field B y 
in the x-direction (except for electrons injected very close to 
the z-axis). This is due to the efficient counterbalance of the 
electric acceleration rate by the Compton energy losses. High 
Compton cooling has a more critical effect on the energy of 
the electrons at the end of the layer. If synchrotron losses are 
negligible (U+ 3> Bq/8tt) and if deviations in the x-direction 
are not too important (xq /L x <C 1), the final energy of the elec- 
tron is given by the balance between the electric force and the 
Compton drag force, i.e. (in the Thomson regime) 



7rad : 



3e/3 lec B ( 



1/2 



(36) 



This expression is the equivalent of the standard radiation re- 
action limited energy derived in Eq. (f3Qb in the case where 
the Compton drag force dominates over the radiation reaction 
force. High inverse Compton cooling severely limits the ac- 
celeration of the particle to extreme energies, but we note that 
this does not affect the concentration of the particle toward 
the reconnection layer. In fact, the Compton drag force helps 
in focusing the orbit toward the xz-plane. 

5. POPULATION STUDIES 

This section is dedicated to the study of dynamics of a pop- 
ulation of particles injected in the reconnection layer. 

5.1. Simulation setup 

We perform an exhaustive exploration of the parameter 
space composed of the initial conditions of the electron, 
namely the Lorentz factor 70, the direction of motion defined 
by the spherical angles 6q = [0,w] (angle with respect to the 
z-axis) and <fio = [0,2ir] (angle with respect to the x-axis in 
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FIG. 10. — An electron trajectory in the yz-plane in a current layer subject to a sinusoidal deformation to mimic the effect of a kink-like instability in the layer, 
with C y = 1 and n z = 5. The location where B K = is shown by the black solid line and the distance ±<5 to this line is indicated by the red dashed lines. Apart 
from B x , all the other parameters of the layer and of the particle are the same as in Fig.[6]but with no guide field. 

the xy-plane), and the coordinates xo = [0,L X ] and yo = [0,5] 
for zo = 0. For each set of initial parameters, the orbit of the 
electron is calcula ted f ollowing the numerical procedure de- 
scribed in Section |4~2l and the value of each of these param- 
eters at the end of the layer is stored in a table and is de- 
noted with a subscript "out". Each particle is followed until it 
crosses the outer boundary of the layer, i.e. \x\ = L x or z = L z 
(see Fig. [12). In practice, yo is set to because the final re- 
sult does not depend much on the value of this parameter (if 

bo I < sy 

From this calculation, we can derive the distributions of the 
accelerated electrons at the end of the reconnection layer, in 
particular the angular and energy distributions. Initially, the 
electrons have isotropic orientations and are uniformly dis- 
tributed in space (along x) and follow a power-law energy 
distribution of index p = 2 such that dNo/d^o oc 7q P with 
7i < 7o < 72- The low-energy cut-off is fixed at 71 = 10 6 and 
corresponds to the Lorentz factor of the bulk particles govern- 
ing the thickness of the layer S = 3.4 x 10 11 cm in Bo = 5 mG 
field strength. The high-energy cut-off is set at 72 = 10 9 , i.e. 
for the particles with a Larmor radius of order the width of the 




FIG. 1 1. — Final Lorentz factor (bottom panel) and deviation angle (top 
panel) of the particle's trajectory as functions of the ratio between the photon 
energy density and the magnetic energy density U±/Ub- Each line corre- 
sponds to different initial distance xq/L^O, 0.1, 0.2, 0.3, 0.5, and 0.9. The 
initial conditions and the properties of the reconnection laye r are the same as 
in Fig. [8] The green dashed line is the solution given in Eq. (36). 



layer L x in the B y 



'■ L x ) = (3 rec Bo magnetic field. In these 



calculations, all large-scale fields are present with a uniform 
E z = /3 rec Z?o, B x following Eq. (TTJ), B y following Eq. © but 
with no magnetic islands or curvature in the layer, and a uni- 
form guide field B z = 0, or 1 times Bo, with the associated 
in-plane electric field components E x and E y . We follow the 
trajectories of 97,200 particles for each simulation. 
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FIG. 12. — This diagram depicts qualitatively the motion of a single elec- 
tron in the layer, injected at some distance Xq from the origin of the axis with a 
Lorentz factor 70. At the end of the layer (z = L z or |.v| = L x ), the Lorentz fac- 
tor of the electron is 7 ou i and its direction of motion is given by the spherical 
angles olll (angle to the z-axis) and <j> onl (angle to the x-axis in the .vy-plane). 



5.2. Outgoing particle angular and energy distributions 

Figure Q~3] shows the outgoing angular distributions of par- 
ticles with respect to 8 oat , dN out /d8 oul summed over cf> out , and 
with respect to <f> out , (W ou t/d0 out summed over 6> out . The dis- 
tribution (W out /d6> out is roughly flat between 6> out = 0° and 90° 
(the distribution exhibits a slight excess between 6* out = 30° 
and 70°) whereas (W out /d0 ou t is concentrated close to <fi out = 
0°, within a few degrees if B z = 0. By symmetry with re- 
spect to the y-axis, if one injects a population of particles at 
xq = [-L x ,0], d/V out /d(9 out is distributed approximatively uni- 
formly between 6> out = -90° and 0°, and (W out /dc/> out is con- 
centrated close to cj) out = 180°. Hence, the reconnection layer 
transforms an initially isotropic distribution of particles into a 
fan beam in the xz-plane. This is associated with the focus- 
ing mechanism due to the alternating reconnecting magnetic 
field and electric field along the z-axis as discussed in detail 
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FIG. 13. — Outgoing angular distributions diVout/dfout (left panels) and dN But /d</> ou t (right panels) of a population of particles (Npat = 97, 200) initially injected 
with a power-law of index —2 with 10 6 < 70 < 10 9 and with an isotropic angular distribution. All particles start at the begining of the layer (y(0) = 0, z(0) = 0) 
and are uniformly spread along the jr-axis between < xq < L x , where 2L X = L z = 4 light-days. The guide field is B z = (top panel), and B z = Bq (bottom panel), 
where Bq = 5 mG. The reconnection rate is ftec =0.1, and Compton losses are neglected. Each line indicates the distribution of particles emitting synchrotron 
photons in the following energy bands: < 10 MeV (green), 10-100 MeV (blue), > 100 MeV (red) and the sum of all bands (dashed line), assuming that electrons 
are bathed in the magnetic field B = Bq. 
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FIG. 14. — Effect of the guide field strength on the outgoing distribution of 
particles emitting synchrotron radiation in the energy bands < 10 MeV (green 
solid line), 10-100 MeV (blue dashed line) and >100 MeV (red dotted- 
dashed line), for electrons bathed in the magnetic field B = Bq. 

in the previous sections. A more precise analysis indicates 
that about 90% of the electrons, regardless of their energy, 
are contained within a solid angle tl out ~ 0.1 steradian, i.e. 
about 1% of the full sphere. A flat distribution of the out- 
going particles with respect to 6> out , AN out /d0 oul as constant, 
implies that the number of particles per unit of solid angle 
diVout/ sin O utdf?out oc 6>~j t is highly concentrated along the z- 
axis, and is the relevant distribution to determine the observ- 
able emission. 

In Figure Q~3] we define three distinct populations of par- 



ticles based on the characteristic synchrotron photon energy 
£syn = (3/2)7Q Ut /jeBo/'« e c that the particles would emit in the 
field Bq once they leave the layer: e syn < 10 MeV (7 out < 3.4 x 
10 8 ), 10 MeV < e syn < 100 MeV (3.4 x 10 8 < 7 out < 10 9 ), and 
e syn > 100 MeV (7 out > 10 9 ). We find that the majority of the 
electrons focused along the z-axis are also the most energetic 
ones. The low-energy particles (7 0Ut < 10 9 ) are dominant at 
high deflection angles 6 out (6> out > 45° for B z = 0). 

A strong guide field decreases the number of particles ac- 
celerated to the highest energies and broadens the angular 
spread of the particles at the end of the layer. Figure [14] 
gives the total number of particles in each energy band as a 
function of the guide field strength. For small guide fields 
(B z /Bo < 0.5), the reconnection layer accelerates up to 50% 
of the particles to the highest energies. Even for a moderately 
high guide field B z /Bq = 0.65, the most energetic particles still 
represent 25% of the total. However, if B z is too strong, com- 
parable to or greater than Bo, then only a few particles can be 
accelerated inside the reconnectio n lay er to the highest ener- 
gies (see the discussion in Section l4~3l ). Most of the particles 
in this case spend most of their time outside the layer, and 
their energy is limited by the radiation reaction limit imposed 
by the external fields (i.e. e syn < 100 MeV). 

The upper panel in Figure Q3] shows the overall energy dis- 
tribution 7o Ut dV ut/d7out of all the electrons at the end of the 
layer for B z /Bq = 0, and 1 . If there is no guide field, the spec- 
trum is well fitted by a single power-law of index +0.5 (this 
result is robust against changes in the initial injection index 
p = 1, 3). The distribution extends almost up to the maximum 
energy that an electron can gain by pure electric linear acceler- 
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FIG. 15. — Outgoing energy distribution of the same population of electr ons 
(top panel) accelerated in the same reconnection layer as described in Fig. 1131 
for B z = (blue solid line), and B z = Bq (blue dashed line). A power-law fit 
to the electron spectrum (with no guide field) dN / d'ye oc is overlaid (ma- 
genta solid line). The maximum energy obtained by pure linear acceleration 
7max = eEoLz I m e c 2 = 3.5 X 10' is shown by the red vertical dotted line. The 
green vertical dotted line indicates the radiation reaction limited energy 7 raI j 
in Bo- The resulting synchrotron emission spectra at the end of the layer are 
also shown (bottom panel) assuming that the particles evolve in a disordered 
magnetic field of strength B = Bq. 

ation 7 max = eEoL z /m e c 2 = 3.5 x 10 9 . The distribution is quite 
narrow, almost monoenergetic. The highest energy particles 
at the end of the layer are mostly composed of the particles 
injected with a moderately high Lorentz factor 70 ~ 10 6 -10 7 , 
because they have enough time to go deep inside the layer 
and they suffer little radiative losses initially. For high guide 
fields, the distribution differs from a single power-law because 
many particles pile-up at the radiation reaction-limited energy 
7 £ w4x 10 s (see Eq. l30l. as mentioned previously. 

5.3. Resulting radiation spectrum 

The lower panel in Figure [15] shows the resulting syn- 
chrotron spectrum emitted by all the particles. For a small 
guide field, it is close to the spectrum of a single elec- 
tron with 7 e = 7 max ps 3.5 x 10 9 , i.e. with frequency de- 
pendence of the flux F v oc v 1 ^ and an exponential cut-off 
at the critical synchrotron frequency u c = 3eB ±~f 2 /47rm e c 
dBlumenthal & Gouldll970l) . For a strong guide field B z = B , 
the radiation spectrum split into two distinct synchrotron com- 
ponents. The low energy component peaks at around 10 MeV 
and corresponds to the synchrotron emission from the parti- 
cles in the radiation reaction limit in Bo = 5 mG field (7 rac j « 
4 x 10 8 ) and /3 lec = 0.1. The high-energy component is fainter 
than the low energy one, peaks at about 1 GeV, and is emitted 
by the particles linearly accelerated to 7 max deep inside the 



layer. 

The synchrotron spectrum is calculated assuming that the 
particles have escaped from the reconnection layer and en- 
tered a region of disordered magnetic field of strength Bq = 
5 mG (the spectrum is averaged over isotropic pitch angles). 
This spectrum represents the instantaneous emission right af- 
ter the end of the layer, assuming that the layer ends abruptljQ. 
For simplicity, this calculation does not take into account the 
photons emitted by the particles during their acceleration in 
the reconnection layer. However, this contribution should be 
small because most of the particles reach the maximum en- 
ergy available in the electrostatic potential, indicating that the 
radiative losses are negligible during the acceleration. This 
statement is not true for a strong guide field (B z > Bo), be- 
cause the particles are stuck in the radiation reaction limit 
during most of their travel along the layer (Fig. [131 . In this 
case, the emission during the acceleration could be as impor- 
tant as or even greater than the radiation emitted at the end 
of the layer. The latter situation should not be considered in 
the context of the Crab Nebula flares, because the maximum 
energy of the photons emitted would not exceed ~ 160 MeV 
contrary to what is observed. The spectrum emitted by a par- 
ticle in the layer does not differ from synchrotron because 
the variation of the Lorentz factor of the particle is small 
over the formation length of a synchrotron photon, which 
is roughly given by the Larmor radius p e divided by 7 e or 
po = m e c 2 /eBo ~ 3.4 x 10 5 cm. (Typically the particle gains a 
small am ount of energy within one half Speiser cycle z\ ^ po, 
see Sect.HTTI) 

Another important issue concerns the calculation of the ra- 
diation spectrum seen by an observer located in a given di- 
rection with respect to the reconnection layer. In the gen- 
eral case, one would need to know the precise structure of the 
fields outside the layer. As the particles leave the layer, they 
are no longer accelerated and their direction of motion may 
isotropize quickly in a B± ~ Bq external magnetic field, hence 
effectively reducing the anisotropy of the emission. How- 
ever, this is not true for the high-energy synchrotron radiation 
above 100 MeV because the emitting particles have energy 
above the radiation reaction limit and cool over a timescale 
shorter than their cyclotron period, i.e. t syn (u>o / 7 e ) < 1, simi- 
larly if /3re C (7 ia d/7e) 2 < 1 (where 7 rad is given in Eq.[30]l. As a 
result, the photons' direction of propagation above 100 MeV 
should be close to the electrons' direction of motion at the end 
of the layer. We use this simplifying assumption in Section[6] 

5.4. Impacts from other parameters 

We ran the same simulations for different thicknesses of 
the reconnection layer ranging from (5 = 3.4xl0 10 cm to 
S = 3.4 x 10 14 cm (« 0.1L X ), but we do not find any signif- 
icant changes in the results. A thick layer tends to acceler- 
ate slightly more particles to the highest energy, especially 
for those that are injected with a high Lorentz factor because 
their Larmor radius becomes comparable to or smaller than 
the thickness 5. These particles then have enough time to fo- 
cus inside the layer where the magnetic field is small. We 
also investigated the effect of inverse Compton energy losses. 
The overall shape of the escaping angular distribution does 
not change, but we note an increase in the den sity of electrons 
at higher angle f? out , as expected (see Section l4~6*l l. The spec- 

1 In reality, the layer may terminate gradually over a distance comparable 
to the length of the layer L z , and the particle would feel a perpendicular field 
strength smaller than Bq. 
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trum remains almost monoenergetic but is notably shifted to 
lower energies if > 0.1 Ub, where the maximum Lorentz 
factor is given by Eq. ( l36*l l. We carried out population studies 
in the presence of magnetic islands in the reconnection layer 
for rii s = 10, 100 and 1000 islands for different amplitudes 

fiy,max/Sy,max = 0.1, 1, and 10. We found that the particles 
cluster quasi-uniformly around the O-points of the islands. 
The outgoing angular distribution of particles is even more 
concentrated at small angle 9 out than without islands. The par- 
ticle energy distribution is also more concentrated at the high- 
est energy j m . dx = eEoL z /m e c 2 = 3.5 x 10 9 because the parti- 
cles are channeled in the z-direction. These features are more 
marked for fi y , m ax/fiy,max ^ 1. Islands increase the anisotropy 
of the particle distribution and the number of extremely high- 
energy particles, but the overall distributions are not changed 
dramatically compared with the case where there are no is- 
lands. Note that these results may not be realistic and follow 
from the assumption of a uniform E z . In a highly dynamic 
plasmoid reconnection regime, E z is expected to be inhomo- 
geneous. 

6. APPLICATION: GAMMA-RAY FLARES IN THE CRAB NEBULA 
6.1. Context and assumptions 

The Crab Nebula is visible throughout the electromagnetic 
spectrum from radio to very-high energy ga mma rays, mostly 
in the form of non-thermal radiation (see lAharonian et al.1 
120041 for a recent compilation of multiwavelength data). The 
broad band Spectral Energy Distribution (SED) is essentially 
composed of a synchrotron bump extending from radio to 
about 100 MeV energy gamma rays and a second bump at- 
tributed mostly to inverse Compton emission emerging above 
100 MeV up to about 100 TeV. The overall SED above UV 
energies is well reproduced by a single radiating popula- 
tion of ultrarelativistic electron-positron pairs injected with 
a power-law energy distribution i mmersed in an average ~ 
100-200 uG magnetic field (e.g.. iHoms & Aharonianl 120041: 
lAbdo et al.H2010t iMever et aLll2010l) . 

The discovery of gamma-ray flares above 100 MeV by 
the gamma-ray space telescopes AGILE and Fermi changed 
this simple steady-state picture. If the variable emission 
above 100 MeV is synchrotron radiation, as it is gener- 
ally believed, and if one associates the duration of the 
flares (?a ~ 1 day) with the synchrotron cooling timescale of 
the emitting particles (t s = 3m e c/2r 2 %B 2 L ), then the gamma 
rays are produced by PeV pairs in a milli-Gauss mag- 
netic field, i.e. in a much stronger field than expected 
from the modeling of the SED, B± 3> 200/iG. This im- 
plies that the particles have to be accelerated beyond the 
radiation reaction limit energy, which is difficult for clas- 
sical models of particle acceleration unless the e mission is 
substantially Doppler boosted toward the observer (|Lvutikovl 
20 lOt iKomissaroy & Lvutikovll201 U feednarek & Idedl201 ll 
Yuan et al.ll2.Ql It iLyutikov et al 1 1201 J). However, optical and 
X-ray observations show that the flow is at mos t mildly rel- 
ativis tic in the nebula (Vfl w/c ~ 0.5, see e.g. iHester et al.l 
120021) . Alternatively, we have shown in Section [4] that 
particles can be efficiently accelerated inside a few light- 
day-long reconnection layer and a milli-Gauss field to ex- 
tremely high energies by the large-scale reconnection elec- 
tric field (Fig. fTsT t. This scenario solves the problem of syn- 
chrotron photons above the classical limit e syn > e* dKirkl 
I2004t lUzdenskv et al.l 1201 lb . In addition, the particles an- 
gular distribution is expected to be extremely beamed along 



the current layer, in particular for the most energetic particles 
emitting above 100 MeV (Fig. [T3l. This could explain why 
about 1% of the nebula radiates 5-10 times more flux above 
100 MeV than the rest of the system. 

Another key feature of the model is the quasi- 
monoenergetic spectrum of the particles at the end of 
the reconnection layer (for B z /Bq < 1, see Fig. [T3J. This 
distribution represents the injection spectrum of fresh parti- 
cles into the nebula and might not always coincide with the 
spectrum of electrons inferred from the observed gamma-ray 
emission averaged over the total duration of the flare. In 
fact, radiative cooling and the (unknown) time-dependent 
injection of these PeV particles can produce broader distribu- 
tions like power-la ws above 100 MeV as reported by Ferm^\ 
dAbdo et al.H2011l) . Reality might be a complex interplay of 
these two effects. The observation of the detailed spectral 
evolution during a flare should provide a powerful diagnostic 
of the injection of PeV particles in the nebula. In addition, 
even a slight change in the orientation of the reconnection 
layer with respect to the observer during the flare, for instance 
due to instabilities in the layer such as tea ring and/or kink 
modes (lBegelmanlll998t iMizuno et~aT1l201 lh . can modify the 
observed emission and explain the short intra-flare variability 
as well as the symmetric shape of t he light-curve profile 
dBalbo et alJ20TltlBuehler eilll20Tlh . 

6.2. Spectral modeling 

For illustrative purposes, we use a simple toy model to re- 
produce qualitatively the gamma-ray emission observed dur- 
ing the September 2010 flare (Figure [T6li. First, the "quies- 
cent" emission of the nebula is obtained from the synchrotron 
emission of a cooled and isotropic distribution of electron- 
positron pairs steadily injected with a power-law energy dis- 
tribution of index 2.2 for 7 e > 10 3 , with an exponential cut-off 
at 7cut = 2 x 10 9 , in a 200 pG magnetic field. On top of the 
quiescent emission, we start injecting, at the beginning of the 
flare (marked as t = 0), another population of particles accel- 
erated in a 4 li ght-day-long reconnection layer. As discussed 
in Section 15.31 it is a good approximation to assume that the 
> 100 MeV synchrotron photons are emitted in the same di- 
rection as the radiating electron's direction of motion at the 
end of the layer. Hence, the observed emission during the flare 
depends strongly on the viewing angle # bs of the reconnec- 
tion layer because of the anisotropy in the outgoing particle 
distribution (Fig.QjJi. The emission is maximum if the layer 
is pointing at the observer (# b s « 0). 

For the spectral modeling, we assume that B z = and 
$obs ~ 0. We extract from Figure Qj] the spectrum of the 
particles contained in the first bin of the angular distribu- 
tion (^out = [0,2.5°]) summed over <f> out . After this opera- 
tion, the spectrum of the particles is composed only of the 
most energetic particles with 10 9 < 7 e < 3 x 10 9 (See Fig.[T6l 
bottom panel). Once injected in the nebula, particles are 
no longer accelerated and are subject only to synchrotron 
cooling in a Bq = 5 mG disordered magnetic field. We ar- 
bitrarily choose that the number of fresh particles supplied 
per unit of time by the layer decreases exponentially with 
time (Wi n j/d7e = dNo/d-f e exp(— f/r), with a decay timescale 
t = 1 day. The particle spectrum time evolution is obtained by 

2 The spectrum of the September 2010 flare reported by the AGILE team 
(averaged over 2 days) is a bit harder than the Fermi spectrum (averaged over 
4 days). This diffe rence could be evide nce of synchrotron cooling of the 
particle distribution ( Vittorini et al. 201 1). 



14 



Cerutti, Uzdensky, & Begelman 




10" 2 10" 1 10° 10 1 10 2 10 3 10" 10 5 
6, [MeV] 



Time evolution electron distribution 




10 B 10 9 10 10 

7e 



FIG. 16. — Top panel: Spectral modeling of the Crab Nebula gamma-ray 
emission during the September 2010 flare. The black dashed line repre- 
sents the "quiescent" spectrum (steady state emission, synchrotron compo- 
nent only) and the red solid line is the synchrotron spectrum of the particles 
accelerated in the reconnection layer averaged over the duration of the flare 
(four days). The magenta dotted lines show the evolution of the flaring spec- 
trum sampled every 12 hours (top line t = and bottom line ( = 3.5 days, see 
the text for more details about the modeling). The black solid line is the sum 
of the qu iescent and the flaring components. SWI FT and Fermi data during 
the flare (Sakamoto et al. 2010; Abdo et al. 201 1) are overplotted with blue 
bowties. Archive data from COMPTEL (Kuiper et al. 2001) and the quies- 
cent Fermi data (only the synchrotron component, Abdo et al. 201 1) are plot- 
ted for comparison. The distance of the Crab Nebula is fixed at 2 kpc. The 
flaring inverse Compton emission is negligible. Bottom panel: Time evolu- 
tion of the radiating electron distribution sampled every 12 hours injected in 
the Crab Nebula during the September 2010 gamma-ray flare. 

solving the Fokker -Planck equation using th e implicit Chang- 
Cooper algorithm (Ch ang & Cooperlll970l) and is shown to- 
gether with the resulting synchrotron spectrum time evolution 
in Figure [16] The electron and gamma-ray spectra soften and 
broaden slightly with time. Nonetheless, little flux is expected 
below 100 MeV compared with the emission from the rest 
of the nebula. Averaged over the duration of the flare, the 
computed spectrum reproduces qualitatively the observations 
(Fig. [16] top panel). We also computed the in verse Comp- 
ton em ission associated with the flare. Following lMever et al.l 
(1201 Oh . the target photon density is composed of the syn- 
chrotron photons from the nebula, thermal emission at 93 K 
from the dust and the 2.7 K cosmic microwave background. 
We find that the expected inverse Compton emission is far too 
weak (by several orders of magnitude) compared to the quies- 
cent Compton emission to increase the very-high energy flux 



noticeably, because the energy density in the background radi- 
ation field is much smaller than the magnetic energy density 
in the Crab Nebula. The inverse Compton drag force has a 
negligible effect on the motion of the particles. In this model, 
the gamma-ray flares are confined to the high-energy band 
> 100 MeV up to a few GeV and do not have any detectable 
counterparts at other wavelengths. 

6.3. Energetics of the flare 

The total number of pairs required in the spectral model- 
ing to account for the September gamma-ray flare luminosity 
above 100 MeV is = 8.5 x 10 37 pairs (if the Crab Neb- 
ula is located at 2 kpc). This number is not corrected for 
the extreme anisotropy of the radiating particle distribution. 
Taking the beaming effect into account, the actual number 
of pairs needed is N ln j = A^°(fV47r) ~ 7 x 10 35 pairs (with 
O « 0.1, see Section [5]). The mean energy of the parti- 
cles injected by the layer into the nebula is (7 e ) « 2.3 x 10 9 . 
Then, the total energy required in these PeV particles is about 
^Fiare ~ 10 39 erg. This energy represents only about 0.2% 
(about 35% if there is no beaming) of the total magnetic en- 
ergy available in the flare region given by 

c 

£mag = -t—E z Bq X 2 (2L x L z )f F i are 
47T 

= 4/3 rec (||) (cfFiare) 3 « 4.4 x 10 41 erg, (37) 

i.e., the Poynting flux passing through the surface of the layer 
(2L X L Z ) = (cfpiare) 2 in the ±y-direction times the duration of 
the flare. 

We conclude that this scenario provides a plausible solution 
to the Crab gamma-ray flares paradox. 

6.4. Possible locations of the reconnection layer 

Our model does not specify the location of the reconnec- 
tion layer responsible for the gamma-ray flares, but one can 
envisage at least three plausible sites: within the pulsar wind 
upstream of the wind termination shock, just outside the ter- 
mination shock in the equatorial region, or outside the ter- 
mination shock in the polar regions. In the case of a mis- 
aligned rotator, the magnetic field in the wind is expected 
to have a toroidal striped-like structure of alternating polarity 
around the equa torial plane (the so-called "striped wind", see 
iKirk et al] 120091 for a review). Magnetic reconnection could 
occur in this configuration and convert the magnetic e nergy 
into particle kinetic energy (ICoronitil 1 990t iMichell 1 9941) . The 
length of the layer would be of orde r the distance to the pulsar. 
However. lLvubarskv & Kirkl (|2001) found that the wind in the 
Crab should reach the termination shock before a substantial 
fraction of the magnetic energy is dissipated, unless the w ind 
contains a high density of pairs (IKirk & Skja eraasen 20031) . 

If the wind remains highly magnetized up to the ter- 
mination shock, the striped structure could be rapidly 
dissipated by shock-driven magnetic reconnect ion in the 
equatorial regions down stream of the shock dLvubarskyl 
2001 iPetri & Lvubarskvl J2007T: iLyubarsky & Livertsl 120081 
Sironi & Spitkovskyil2011l) . similarly to what may happen in 
the nonrelativistic case of th e solar system clo se to the he- 
liospheric termination shock (Drake et al. 2010). The current 
layer could be of order the termination shock radius, ~ 0. 1 pc 
in the Crab, which is too long to explain the short variability 
timescale of the flares, but fragmentation of the current sheet 
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could lead to more rapid variability. The amount of magnetic 
energy available for dissipation in this case would be limited 
by the pulsar misalignment angle, with a high efficiency pos- 
sible only if the angle is close to 90°. 

The third possibility is that rapid reconnection is triggered 
by current-driven MHD instabilities downstream of the shock. 
If the magnetic field is dynamically important, the postshock 
wind will try to relax to a configuration in which gas pres- 
sure forces balance the stresses due to the predominently 
toroidal magnetic field (IBegelman & Lilll992l) . Such a sys- 
tem is highly unstable to kink instabilities (and slightly less 
unstable to pinch instabilities) in regions where the plasma 
beta paramete r is > 0(1), on time scales of order the Alfven 
crossing time jBegelmarJ1998tlMizuno et al.l201 ll) . The non- 
linear development of the instability leads to a complex sys- 
tem of current sheets (O'Neill et al., in preparation) that could 
rapidly dissipate most of the magnetic energy at mid- to polar 
latitudes. The dissipation of the magnetic energy in the polar 
region could power the X-ray jet in the Crab. New observa- 
tions with Chandra show that the southern jet has moved since 
2000 dWeisskopfll2011l) . perhaps related to the kink instabil- 
ity. In addition to the formation of current layers in the polar 
regions, it is expected that the magnetic field strength is am- 
plified close to the rotational axis by the z-pinch, as required 
to power the flares. 

In addition to setting up the conditions for extreme parti- 
cle acceleration, the establishment of reconnection sites due 
to current-driven instabilities could resolve the long-standing 
"sigma problem". Hot gas in polar regions where the mag- 
netic field has dissipated could exhaust to larger radii, allow- 
ing more highly magnetized gas from equatorial latitudes to 
fill the void. As it does so, this fresh gas will also go unsta- 
ble, leading to the destruction of its magnetic energy, and so 
forth. Thus, current-driven instabilities at high latitudes could 
set up a conveyor belt leading to the dissipation of most of 
the magnetic energy introduced by the shocked pulsar wind. 
Such a process could reconcile the relatively high magnetic 
field (and inferred wind a ~ O(l)) required by the particle 
acceler ation model with the low a required by mo dels of the 
nebula dRees & Gunnll 1974 IBegelman & Lilll992l) . 

7. CONCLUSION 

The existence of synchrotron photons above 100 MeV emit- 
ted by <? + /e~ pairs requires the non-standard condition that the 
electric field exceeds the magnetic field in the acceleration re- 
gion. This condition can be fulfilled in regions where ideal 
MHD approximation breaks down, such as reconnection sites. 

Relativistic test-particle simulations presented in this pa- 
per show that high-energy particles moving across a thin re- 
connection layer are accelerated by the electric field and fo- 
cus rapidly deep into the layer where E > B. This mecha- 
nism transforms any initial distribution of particles into an 
extremely beamed distribution in the reconnection layer. In 
addition, our population studies show that the energy distri- 
bution of the particles at the end of the layer piles up to the 
maximum potential drop energy available, eEL z , suggesting 
that the reconnection layer acts almost as a pure linear ac- 
celerator. Because the perpendicular magnetic field is small 
deep inside the layer, individual particles can exceed the nom- 
inal radiation reaction-limited energy and emit > 100 MeV 
synchrotron photons provided that the layer is long enough. 
We note that a strong guide field and the associated in-plane 
electric field tend to defocus the particles outside the layer, 
decreasing the fraction of particles accelerated to the high- 



est energies. We find also that static magnetic islands in the 
layer, generated for instance by tearing instabilities, tends to 
trap and focus the particles toward the O-points. If the layer 
is bent by i.e. kink instabilities, particles' trajectories follow 
the curvature of the layer. The final energy of the particles 
and the general focusing mechanism described in this article 
remain unchanged by these deformations as long as the curva- 
ture length-scale is not too short compared with the distance 
traveled by the particles between two crossings through the 
layer. 

These key features of the model (i.e. generation of a fan 
beam of particles with energy <~ eEL z ) are essential to explain 
the apparent high number of PeV electrons required to power 
the Crab Nebula gamma-ray flares. For illustrative purposes, 
we applied the model to the September 2010 flare. Observa- 
tions can be qualitatively reproduced assuming a few light- 
day-long reconnection layer and a few milli-Gauss reconnect- 
ing magnetic field. The resulting spectrum is close to the syn- 
chrotron spectrum of a single electron. Cooling and time- 
dependent injection of fresh particles in the nebula during the 
flare can broaden and soften the spectrum averaged over the 
duration of the flare. The strong beaming of the electron dis- 
tribution emitting the high-energy radiation accounts for the 
overall luminosity of the flare, if the reconnection layer points 
toward the observer. This implies that gamma-ray flares could 
happen frequently in the nebula, but only those occasionally 
pointing at us are detected. The inverse Compton emission 
associated with the flare is far too low to be observed and no 
detectable counterparts of the flare at other wavelengths are 
expected in this model. 

The mechanism of particle acceleration described in this 
article may be at work in other astrophysical objects, for in- 
stance in active galactic nuclei (AGN) jets (Nalewajko et al., 
in preparation). A few blazars flaring in TeV gamma rays 
exhibit variability timescales (~ minutes) much shorter than 
the light-crossing time of the Schwarzs child radius of the cen- 
tral supermassive black hole (~ hour) dAharonian et aLll2007t 
I Albert et all 120071: lAleksic et all 120 111) , suggesting that par- 
ticles are quickly and efficiently accelerated in a small vol- 
ume of the relativistic jet. The acceleration and, more im- 
portantly, the focusing of the particles within a reconnection 
layer could account for strong beaming of the TeV flare emis- 
sion. In the environment of AGN jets, the effect of the inverse 
Compton drag force becomes important and would decrease 
the maximum energy reached by the particles at the end of the 
layer. Nevertheless, we found that the distribution of particles 
should remain extremely anisotr opic. Note that this model 
differs from the "mini- jet" model dGiannios et al .1120091 120 10t 
iNalewaiko et al.ll2.01 ll) in which particles are assumed to be 
beamed just by the relativistic bulk motion of the reconnec- 
tion outflow (moving perpendicular to the electric field, in the 
jet comoving frame) expected in a Poynting flux-dominated 
jet (Tbuik ~ er 1//2 , where a 3> 1 is the ratio of the Poynting 
flux to the particle kinetic energy flux). In the Crab Nebula, 
mini-jets are unlikely to form because a is at most of order 
unity. 

The test-particle simulations presented in this work are not 
fully satisfactory, because many uncertainties remain regard- 
ing the structure of the electromagnetic fields inside and out- 
side the reconnection layer. Interactions between particles are 
also neglected in this study. It is unclear what the effect of 
extremely high-energy particles on the dynamics of the re- 
connection layer would be. A self-consistent and more realis- 
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tic approach to this problem would be to use PIC simulations 
with ultrarelativistic electron-positron pairs, including the ra- 
diation reaction force as pioneered by iJaroschek & Hoshinol 
(120091) . Another challenge is to calculate accurately the radi- 
ation emitted by the particles during the acceleration process 
inside and outside the reconnection layer (time dependent an- 
gular and energy distributions). These issues are left for future 



investigations. 
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